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Abstract 

Terrestrial Gamma-Ray Flashes (TGFs) are very short bursts of high energy photons and elec- 
trons originating in Earth's atmosphere. We present here a localization study of TGFs carried 
out at gamma-ray energies above 20 MeV based on an innovative event selection method. We use 
the AGILE satellite Silicon Tracker data that for the first time have been correlated with TGFs 
detected by the AGILE Mini-Calorimeter. We detect 8 TGFs with gamma-ray photons of ener- 
gies above 20 MeV localized by the AGILE gamma-ray imager with an accuracy of ~ 5 — 10° at 
50 MeV. Remarkably, all TGF-associated gamma rays are compatible with a terrestrial production 
site closer to the sub-satellite point than 400 km. Considering that our gamma rays reach the AG- 
ILE satellite at 540 km altitude with limited scattering or attenuation, our measurements provide 
the first precise direct localization of TGFs from space. 
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INTRODUCTION 



Earth's atmospheric events associated with strong thunderstorms have been in recent 
years observed to be the site of very efficient particle acceleration and gamma-ray emis- 
sion at MeV energies and above jl-4|. Of particular interest are the so-called Terrestrial 
Gamma- Ray Flashes (TGFs) that are very short (lasting up to a few milliseconds) bursts of 
high-energy photons above 100 keV, first detected by the BATSE instrument on board the 
Compton Observatory TGFs have been associated with strong thunderstorms mostly 
concentrated in the Earth's equatorial and tropical regions {2, ^J. TGFs are widely believed 
to be produced by Bremsstrahlung in the atmospheric layers by a population of runaway 
electrons accelerated to relativistic energies by strong electric fields inside or above thun- 
derclouds. The secondaries generated during the acceleration process can be accelerated 
as well driving an avalanche multiplication [6j, commonly referred to as Relativistic Run- 
away Electron Avalanche (RREA). However, the RREA mechanism alone is not sufficient 
to explain the rich phenomenology of TGFs, especially the observed fluence, and there is no 
consensus yet on the underlying physical conditions, production sites, radiation efficiencies 
and maximal energies. An interesting possibility to overcome some of these difficulties is 
the relativistic feedback mechanism [7J, |8j , which predicts an avalanche multiplication factor 
and characteristic discharge time compatible with the observed fluence and time profile of 
TGFs. Concerning the maximal energy, the original BATSE detection of TGFs up to a few 
MeV Q, EJ was superseded by the RHESSI detection up to 20 MeV Q, Q. Recently, the 
AGILE satellite showed that TGF spectrum extends well above 20 MeV (Tavani et al., 
2010, submitted to Nature), as confirmed also by the Fermz-GBM detector 12]. 



AGILE 13| is a mission of the Italian Space Agency (ASI) dedicated to astrophysics in 
the gamma-ray energy range 30 MeV - 30 GeV, with a monitor in the X-ray band 18 keV - 
60 keV [ijj], operating since April 2007 in a low inclination (2.5°) Low-Earth Orbit at 540 km 
altitude. The AGILE Gamma-Ray Imaging Detector (GRID) is a pair-tracking telescope 
based on a tungsten-silicon tracker [15j. The imaging principle is based on the reconstruction 
of the tracks left in the silicon detection planes by the electron-positron pairs produced by 
the primary photon converting mainly in the tracker tungsten planes. A Mini-Calorimeter 



(MCAL) 16j, based on CsI(Tl) scintillating bars for the detection o: 



range 300 keV - 100 MeV, and a plastic anti-coincidence detector 



gamma-rays in the 



17] complete the high- 



energy instrument. MCAL can work also as an independent gamma-ray transient detector 
with a dedicated trigger logic acting on several time scales spanning four orders of magnitude 



between 290 fis and 8 seconds 18|, ll9j . Thanks to its flexible trigger logic on sub- millisecond 



time scales, MCAL proved to be a very efficient instrument for TGF detection. The average 
MCAL detection rate is ~ f0 TGFs/month, with the current selection criteria 11]. 

Up to now TGF observations have only been reported by space instruments with no or 
quite limited on-board imaging capabilities (e.g. BATSE and Fermi-GBM). The aim of 
this Letter is to provide a first accurate localization of TGFs from a space instrument and 
to study their very significant high-energy tail of emission above 20 MeV. 



AGILE-GRID DETECTION OF TGFS 

In the period between June 2008 and December 2009 the MCAL instrument triggered 



119 bursts identified as TGFs according to the selection criteria discussed in [ll|. For each 
of these bursts, the GRID dataset was searched for quasi-simultaneous gamma-ray events 
within a 200 ms time- window centered at the TGF start time T , defined as the time of the 
first MCAL-photon associated with the TGF. Figure [T] shows the cumulative histogram of 
the arrival times of the GRID events with respect to T obtained by summing all the 119 
bursts. A peak in the distribution is evident for the 2 ms time bin immediately following To. 
This peak includes 13 events, and the probability for it to be a statistical fluctuation (13 
events or higher) is 6.5T0 -10 if we assume that GRID events are not correlated to TGFs and 
are distributed according to the Poisson law with the measured average rate of 5.1 counts/s. 
All these GRID events take place during the TGF emission time interval estimated from 
MCAL data only. 

A standard direction reconstruction of every GRID event is performed on board by means 



of a GRID-adapted Kalman filter [20j, |21[ : the event is discarded if the incoming direction 
is found to be within 70° from the Earth center. This procedure, called albedo filtering, is 
aimed at the rejection of the Earth gamma-ray albedo photons which are one of the most 
significant contribution to the gamma-ray background. This procedure applies to all GRID 
events considered in our analysis, except for the events obtained during a period of about 
100 days (1495 orbits) during which the albedo filtering was disabled for test purposes. For 
any given track in the GRID, the albedo filtering assumes an incoming direction compatible 
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FIG. 1. Histogram of the arrival times with respect to the TGF MCAL-start time To for the 
gamma-rays detected by the AGILE gamma-ray imaging detector (GRID). 



TABLE I. Selected events properties. 
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a Start time of the TGF, defined as the time of arrival of the first photon 
b TGF duration 

c Time difference between GRID event time and To 

d F: the event direction intersects the Earth. R: the reversed event direction intersects the Earth. 
c Angle between event direction and the Nadir. A 5.8° uncertainty can be considered for an average 



60 MeV photon energy. 

f Distance of the event direction projected to the Earth surface from the satellite footprint. The 

uncertainty due to the error on direction reconstruction is between 50-90 km for an average 60 MeV 

photon. 
s Albedo filtering disabled 

with the satellite field of view and cannot easily discriminate whether a photon producing a 
single or quasi-single (i.e. an event configuration of ambiguous topology) track in the GRID 
originated from the opposite direction. This process is especially important for low energy 
(few tens of MeV) photons which tend to produce single or quasi-single tracks with a few 
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detection planes involved. 



All the selected GRID events were processed by the AGILE standard analysis pipeline 
which iteratively fits the tracks by means of a custom Kalman filter [21] in order to extract 
the most probable energy and direction of the incoming photon. The nominal uncertainty 
on energy estimates is a factor ~ 2, mainly due to the fact that MCAL is only 1.5 radiation 
lengths thick and cannot provide full energy containment for electromagnetic showers. The 
angular resolution (68% containment radius) is 3.5° at 100 MeV [3] and scales approxi- 
mately as the inverse of energy, being dominated by the multiple scattering effects which 
affect mostly the less energetic tracks. Using the satellite position provided by the GPS 
system, it is checked whether the event's incoming direction is compatible with an origin on 
the Earth surface, parameterized as the World Geodetic System WGS84 ellipsoid. If this 
is true, the event is flagged as a Forward (F) event. If no intersection is found, the event 
is processed again with a non-standard algorithm that assumes that the track vertex is in 
the bottom planes of the tracker and the tracks are allowed to develop from the bottom to 
the front planes. This algorithm returns new estimates of energy and incoming direction. 
It is then checked again whether the new incoming direction is compatible with the Earth's 
surface. If this is the case, the event is flagged as a Reverse (R) event. Nine events out 
of 13 are classified as either F or R events. For these 9 events, Table [J reports the main 
properties of the TGFs and associated GRID events. It is worth noting that all the 3 F 
events have been detected during observations with the albedo filtering disabled. The effects 



of a TGF production altitude between 15-40 km 10j, |22j have not been considered at this 
stage because they are much smaller than the uncertainties due to the instrument angular 
resolution. Among the four GRID events with directions not compatible with the Earth's 
surface, three of them come from directions very close to the Earth's limb, which make them 
classifiable as albedo photons passing the on-board albedo filtering step. 

All TGFs in the selected sample have one associated GRID event exept TGF 12809-19, 
for which two closely spaced GRID events were detected. For this TGF, Fig. [2] shows the 
topology of the tracks associated to the two photons detected in the tracker. The tracks 
shown join the hits detected in the silicon tracker planes (the planes are not shown for 
clarity). For both events, a track vertex can be found indicating that the events are coming 
from the bottom of the instrument. 



6 




FIG. 2. 3D view of the tracks reconstructed for the two GRID events detected in association with 
TGF 12809-19. The MCAL detector with the detected events is shown at the bottom of the figure. 
The AGILE pointing direction is orthogonal to the MCAL detection planes, toward the top of the 
figure. 

RESULTS AND DISCUSSION 

The incoming directions of the 9 selected events appear to be clustered close to the sub- 
satellite point, with an average (maximum) angle (the angle between the photon direction 
and the satellite nadir) of 25.4° (35.1°) and distance to the subsatellite point of 260 km 
(390 km). All 9 events are thus contained within a 1.14 sr solid angle, a factor 3.4 smaller 
than the solid angle subtended by the Earth at the satellite altitude of 540 km, which cor- 
responds to a maximum visibility projected distance radius of ~ 2600 km from the satellite 
footprint. 

Figure [3] shows the geographical location of the satellite for the considered TGFs as 
well as the location of the intersection of the tracked events direction with the Earth. It 
is remarkable to note that the two GRID events associated to TGF 12809-19, both with 
energies of about 40 MeV, come from very close directions, separated only by a 5.4° angle, 
within the instrument angular resolution at these energies. The time delay between the two 
events is only 264 fis. Considering that the dead-time for GRID data acquisition is 200 /is, 
it is possible that high energy measurements for this TGF, as well as for others, is hampered 
by the dead-time. Figure H] shows the scatterplot of the GRID events projection with respect 
to the AGILE footprint, and the distribution of the occurrence density vs. distance from 




FIG. 3. Location of the satellite footprints (filled circles) and of the projection of the associated 
tracker events to the Earth ellipsoid (hollow circles). The inset is a zoom over the region of 
TGF 12809-19. The dashed ellipses represent the la error circle for the associated tracker events. 
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FIG. 4. (a) scatterplot of the GRID events projection with respect to the AGILE footprint. The 
dashed circles are 300 and 600 km in radius. The full AGILE visibility region has a 2600 km radius, 
(b) occurrence density vs. distance from footprint (each bin has been divided by the subtended 
area in km 2 ). 

footprint (each bin has been divided by the subtended area in km 2 ). 

Gamma-ray photons can be Compton scattered in their path from the TGF source to the 
satellite. We calculated the scattering probability and an gula r distribution of the resulting 



gamma-ray flux for a TGF source at 15-40 km altitude 
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221 ]. using the photon cross- 



sections provided by the NIST XCOM database |23j and assuming an air density vertical 
profile scaling exponentially with a length scale of 7 km. Above 40 km, the probability for a 
photon of energy larger than 20 MeV to be Compton scattered results lower than 3%. If a 
Compton interaction takes place in the 15-40 km altitude interval, the resulting positional 
smearing will be within the angular resolution of the GRID reconstruction, considering also 
that small scattering angles are favoured at high energy. Furthermore, we also considered 
the possibility of the production of secondaries from below originating either by Compton 



scattering or pair production in the spacecraft or in MCAL, which can contribute to our 
R-detections. In both cases, the resulting angular smearing is of the same order of the GRID 
angular resolution. 

Our detection of TGFs with an imaging gamma-ray detector is important for several 
reasons: (1) it provides the first accurate localization (within a few degrees) of TGFs from 
space; (2) it reinforces the AGILE discovery that a significant high-energy component of 
emission is produced by TGFs well above 20 MeV; (3) it shows that the TGF high-energy 
emission is detected by satellites in LEO orbits from a relatively small region within 300-400 
km from the satellite footprint. The AGILE detection of gamma-ray photons above 50 MeV 
from TGFs (the average energy of the 9 GRID events is 60 MeV) strongly constrains the the- 
oretical models. For example, the detection of a given energy sets a lower limit in the electric 
potential difference involved. In the relativistic feedback scenario 3], where the maximum 
possible static electric field is limited by the avalanche mechanism, this correspond to a lower 
limit in the column depth of the avalanche region hence, given a production altitude, sets 
a lower limit in the extension of the avalanche region. Moreover, high energy electrons are 
expected to be very well aligned with the electric field 24J. Considering the small angular 
scattering of Bremsstrahlung photons for highly relativistic electrons, the incoming direction 
of high energy photons marks well the electric field orientation at the source. In this sense 
our localization technique provides also a diagnostic tool for the electric field at the source 
region. The observed clustering of TGFs close to the subsatellite point is in agreement with 
other independent determinations of TGF locations based on ground measurements (sf er- 
ics) simultaneously obtained with space detections by the RHESSI satellite |25l-l27|. Future 
investigations will determine whether the apparently narrow cone of detected gamma-rays is 
due to beamed emission intrinsic to the TGF source or caused by a selection effect favored 
by absorption and Compton scattering in the atmosphere. 

AGILE is a mission of the Italian Space Agency (ASI), with co-participation of INAF (Is- 
tituto Nazionale di Astrofisica) and INFN (Istituto Nazionale di Fisica Nucleare). Research 
partially funded through the ASI contract n. 1/089/06/2. 
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